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The thermal relaxation time (r Kcc ) for the degenerate electron plasma has been calculated by 
incorporating non-Fermi liquid (NFL) corrections both for the thermal conductivity and specific 
heat capacity. Perturbative results are presented by making expansion in T/itld with next to leading 
order corrections. It is seen that unlike the normal Fermi liquid (FL) result where r Kcc oc l/T 2 , NFL 
corrections in leading order (LO) changes the temperature dependence of r Kee to l/T. Incorporation 
of the phase space correction driven by the medium modified Fermion dispersion relation increases 
^— ^ the relaxation time further. 
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C/5 I I. INTRODUCTION 

CSl . Determination of thermal relaxation time of degenerate electron matter has been a subject of serious investigation 
for the last several decades. Application of such studies encompass broad areas spreading across various disciplines 
like metals, semiconductors, astrophysical objects like white dwarfs or neutron stars to name a few. Our focus here is 
to study the heat conduction in neutron star. In particular, we determine the thermal relaxation time of degenerate 
electron system at high density which is relevant for heat transfer from the crust of a star to the core[l|-[3j]. 

It is known that when a new star is born following a supernova explosion, large amount of neutrinos are emitted 
immediately from the core resulting in colder core and a hotter crust, thus a temperature gradient is set up. Then the 
thermal energy gradually flows inward by heat conduction which alternatively might be viewed as the propagation 
of the cooling waves from the center towards the surface leading to thermalization One of the subjects of 
contemporary research in astrophysics has been the estimation of this thermalization time scale or estimation of the 
thermal relaxation time. The investigation what we pursue here is relevant in the context of neutron star as we know 
that degenerate electron gas and positively charged ion constitute the envelop of the crust of the neutron stars. 
There exists several calculations where heat conduction has been studied extensively. In these investigations major 
1 contributions have been seen to originate from the electron-ion scattering 0, The contribution of the electron- 
electron scattering, in contrast, have been found to be of limited importance [J|. Recently this problem has been 
■ revisited and it was seen that such conclusions are true only when one considers the charge charge interaction and 
, neglects the current-current interaction completely. This is a valid approximation in dealing with the ions, but might 
not be justified for the electron-electron scattering where at high density the magnetic interaction becomes important 
043 ■ Reference @, in particular, deals with the calculation of heat conductivity where it has been shown that at 
high density, due to strong magnetic interaction, the electron-electron collision frequency become larger than the 
1 electron-ion collision frequency reducing the heat conductivity. The other point is to note that almost all these 
H . calculations treat the degenerate electron matter as ideal Fermi liquid and treat the electron-electron and electron-ion 
scattering non-relativistically restricting to the electric sector. But at high density for the electrons with momentum 
close to the Fermi momentum, since relativistic effects become important, the magnetic interaction can no longer be 
neglected. It is now known that with the inclusion of the transverse interaction the normal Fermi liquid description 
breaks down due to the vanishing of the electron propagator near the Fermi surface. This can be attributed to the 
absence of static screening of the magnetic photon [8|. Several investigations have been performed in recent years 
where incorporation of such corrections have been seen to have serious implications on various physical quantities like 
pressure, entropy, viscosity or quantities like drag and diffusion coefficients P4l^|. Non- Fermi liquid behavior for the 
neutrino emissivity or the neutron mean free path have also been studied extensively [l3l - fl5| . In all these calculations 
such NFL corrections have been observed to be significant compared to the Fermi liquid results. 

In this work, therefore we incorporate NFL corrections while estimating the thermal relaxation time (r Ke ) of 
degenerate electron plasma. Such estimation requires knowledge of both the thermal conductivity (ft e ) and specific 
heat (c v ) where this correction has to be included consistently for a given order. Furthermore in dense matter the 
quasiparticlc dispersion characteristics change which modify the density of states too. Inclusion of this, as we shall 
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see, also modifies the results significantly both for n e and r Ke . To derive analytical expressions for these quantities 
however, we make perturbative expansion in T/mu, where, T is the plasma temperature and mjj is the Debye mass. 

The plan of the paper is as follows. In section II we develop the formalism part which incorporates the Boltzmann 
equation, the screening mechanism of long-range interactions and the evaluation of the thermal relaxation time. In 
subsection A the results of leading order thermal conductivity and the thermal relaxation time have been discussed 
and in the subsection B next to leading order NFL correction of the thermal conductivity and the thermal relaxation 
time have been included followed by summary and conclusion. 



II. FORMALISM 



We consider degenerate electron gas where the electrons constitute an almost ideal and uniform Fermi gas and 
collide between themselves. We aim to calculate the thermal relaxation time of degenerate electron gas with the 
medium modified phase space factor. The characteristic relaxation time for thermal conduction t k can be defined as 
follows 

r K = -■ (1) 

c v 

In case of the strongly degenerate electron gas the electron thermal conductivity (« e ) can be expressed as follows, 

Jt 

K e = -= — , v e = v ei + v ee , (2) 
1 v e 

where, Jt is the thermal current and v e is the total effective collision frequency. This one is the sum of the partial 
collision frequencies, i.e, electron-ion (v e i) and electron-electron (v ee ) scattering rate. Evidently n e is related to K ee 
and n e i ta|, 



(3) 
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Hence, derivation of K ee in turn requires the knowledge of v ee . It is clear from the denominator of the above equation 
that the heat conduction becomes difficult when collision frequency increases. 

To evaluate re ee we appeal to the Boltzmann equation which describes the kinetics for the individual fcrmion 
component 



0_ 

at 



Vp.V P + P.Vp / P = -C[/ P ], (4) 



here, p is the momentum of the quasiparticle, F is the external force, v p is the velocity of the heat carrier and f p 
is the distribution function of electrons. The collision-integral on the right-hand-side (RHS) is given by the rate of 
fermions scattering in and out of the state with momentum p by scattering on the other fermions with momentum p . 
In the presence of a weak stationary temperature gradient and absence of any external force the Boltzmann equation 
takes the following form, 

Vp.V r /p = -C[/ p ]. (5) 

Now, due to presence of a weak temperature gradient these Fermi-Dirac distribution functions deviate from equilibrium 
distribution functions /,-, which we write as, 

l = h + ^^ (6) 



where. 



/i = < [exp(^)+l<> . (7) 
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£ is the particle energy, /j, is the chemical potential and T is the temperature. Clearly the second term with $ 
measures the deviation from equilibrium. The collision integral can be written as follows, 



C[f P ] 



Mi ± /;)(i ± f k ) - ± / p )(i ± f k )(2n) 4 s 4 ( P + P - k - k')\Mf 



(8) 



\M\ 2 is the squared matrix element for the scattering process 12 — } 34. The ± sign include stimulated emission and 
Pauli blocking. Using the standard linearization procedure from the Eqs.([5|) and (|6]) we obtain the equation for $, 



/p(Wp)(e P -M>V = ^ / f P f k (l±f p )(l±f k )(2^6\p + k-p -k )\M\\% + $ k -%' -**')• (9) 

J P ',fc,fc' 

The above equation can be written in the form \X) = I\<&), where, \X) = (e p — n)v p and / is the integral operator. 
The thermal conductivity n ee is given by the maximum of the following equation and has already been discussed in 



id, 



(X\$) 2 



(10) 



|-) denotes an inner product, the quantity 



(A'|^) 2 



T(*|/|*) 

way to define K ee in the Eq.([3]). Hence, one can write 



is minimal for ^ = $ with the minimal value k pp . This is another 
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Vzfp(l - f P )% 



; f p f k (l ± f p )(l ± f' k )(2i:) A 5\p + k-p -k )\M\ 



(ii) 



the term in the first bracket in the denominator is the thermal current Jt- $ should be determined by minimizingthe 
Eg. (1101) and the minimal value is 'J. But for the present purpose here we consider the simplest trial function 



ty p tx (e p - /j,)v z 



(12) 



The above trial function can now be inserted in the Eq. ([TTj) and the term in the bracket can be averaged over the z 
axis keeping x and (f> fixed, where, (f> is the azimuthal angle between v p and v k . After averaging we obtain [f|, 



(* p + * k - - * fc ,) 2 = -c 2 (l - x 2 )(l - cos0). 
To proceed further one needs to know the interaction, considering only the electron-electron scattering [f| , 



(13) 
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The above matrix amplitude squared is infrared divergent in the t — > limit. To circumvent this problem one needs to 
consider the medium modified screening effect. We in this paper use the Braaten and Pisarski's prescription [T^| . In 
this method in the lower momentum region hard-dense-loop (HDL) resummed propagator [IUGjl is incorporated and 
in the higher momentum region one can use the bare propagator. Using the in medium HDL resummed propagator 
the matrix amplitude squared for the electron electron interaction takes the following form, 



IMI 2 = 32e 4 



(1 — x 2 )cos(j> 



(q 2 + n L ) (q 2 



n T ) 



(15) 



where, and are the longitudinal and the transverse part of the dressed photon propagator and in the HDL 
framework the structure of these two polarization tensors are the following, 



n L (g, uj) = m 2 D XL, n T (g, w) = m 2 D XT, 



(16) 



where, 



XL 



XT 



1 - -In 
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x- 1 
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In the above expressions x = qo/q and is the Debye mass m 2 D = e 2 fj, 2 /tt 2 . In case of the thermal conductivity 
due to presence of cos0 term in the numerator of the Eq. (|llj) , the cross term of the matrix amplitude squared does 
not vanish, hence, on retaining the cross term the matrix amplitude squared becomes, 



IMI 



32e 4 



2g 4 cos</ 



q 2 cos 2 ' 



(I 2 + ™D? (q 2 + m 2 D ) («« + £*g&) + g« + 



(18) 



Now, we first compute the denominator of the Eq. (|ll[) . The denominator is the thermal current Jt as indicated 
earlier and is given by, 
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( e P - m) 
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v z f P (l - f P )% 



H 2 T 2 



(19) 



To proceed further to evaluate the numerator in the Eq. (|llj) it is convenient to introduce a dummy integration variable 
uj. We write the energy conserving delta function as, 



6{e p + e k - e p - e k ) 



doj5(u) + 



)S(u + e k -e k ). 



(20) 



Evaluating q = p — p in terms of p, q and cos9 pq and defining t = uj 2 — q 2 wc find, 



5(ui + e p — e ) = — S co&0 pq — 
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1kq 
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5{lj + e k - e k ) 
using the above delta functions we obtain, 

|fe J ' dpdk f e J ek (l± /(e p + w ))(l± f^-^dqduj^il- x 2 ) ' 



2n 



;i -cos0)|Af| 5 



(21) 



(22) 



The inclusion of the electron self-energy in the dispersion relation changes the phase space factor, which, in turn, 
changes the energy integral in the above equation. We write the momentum integration in the phase-space factor of 
the Eq.flMI) as, 



dk 

de k 



1 



<9Re£\ 



did J ' 

where, we have assumed that the quasiparticle energy uj obeys the following dispersion relation, 

u = e k — Rc£(u;, k). 



(23) 



(24) 



Only the real part of the self-energy has been taken into consideration as the imaginary part turns out to be negligible 
in comparison with the real part. In the free case dk/de k gives the inverse fermion velocity. 
From the Eg. (|11|> we can now write the expression for as, 
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(25) 



where, 



I Kee (T/m D ) = 



diu 



4T2 



(sinh^) 
1 



dx 



2n 



x 2 (\ - x 2 )(l - cosq 



1 + (xniD/uj) 2 XL(x) 1 + {xrriD I uj) 2 xt{x) / (1 - x 2 ) 



(26) 
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In the Eg. (f2l)]) the major contribution comes from the small angle approximation i.e the small value of x dominates 
(x ~ uj/rriD or equivalently q ~ mo). We have to consider the small x behavior of Hl,t, we approximate the 
polarization function in the small x (\x\ << 1) limit to obtain, 

7TT* 

XL = l + 0(x), XT ( x )= l — +0 (x 2 ). (27) 

In the next two subsections we present the results of the thermal relaxation time both for the leading and higher 
orders. 



A. Leading order thermal relaxation time 



In this section we derive the LO result of the thermal relaxation time in the degenerate electron plasma present in 
the outer crust of the neutron star. It might be recalled that in the FL theory the magnetic contribution is suppressed 
compared to the electric one. In this domain n ee varies inversely with T and c v cx T. The thermal relaxation time 
on the other hand is the ratio of these two quantities. Therefore, in FL theory or in absence of transverse interaction 
t Kcc oc 1/T 2 . Now, we incorporate relativistic effects in it in the small energy transfer region (cj << mjj, T « mjj) 
where transverse, weak dynamical screening effect becomes important. For this, transverse interaction has been 
incorporated in the thermal conductivity and medium modified phase-space factor has also been included. Now, to 
have the LO expression of T Kee we have to evaluate I Kcc (T/mu) in the Eq. (|26p . In this region where uj << tud, 
T < < vein the upper limit of the integral can be sent to infinity. Electric interaction in this region gives the following 
contribution, 

f°° du u 2 [°° , x 2 (1 - x 2 ) 

II = I 9 / dx 

Jo w 4T 2 (sinh^) Jo ( x , ™*nf s ' 



where, A is the small temperature expansion parameter T '/mo- The magnetic interaction on the other hand gives, 

du uj 2 r , x 2 (i - x 2 ) 

dx— 



lo u AT 2 (sinh#) Jo l -p 
~ 2A 2 C(3). (29) 

The higher order contribution ((T/m^)) 10 / 3 ) in the above equation can be neglected. For the cross term we find that, 

du uj 2 r , x 2 (i - x 2 ) 

Il,t = I n I dx- 



3 " v 3 / V 3 

From the above Eqs. ([2^1 [2^1 [301) we see that the term coming from the magnetic sector dominates over the electric 
and the crossed terms. After the angular integration we now focus on the momentum integral. The momentum 
integration gives, 

/ dpdkf e J Ck (l - f {ep+u) )(l - f (ek - u) ) = (1 + 7 ) T 2 V* (31) 

where, we have used Eqs. ((23l [24]) . The details of this integration is given in the appendix. In the above integral 7 
represents the medium effect in the phase space factor. In the LO, 7 receives the logarithmic correction from the 
quasiparticle self-energy, 

'Ay/2: 



S(e,*) = e' m {-V[togF^) + ll}, (32) 
1 127H?n L V 7re / J J 
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where, e is chosen to be (e^ — /x) and in the low temperature limit (e^ — /i) ~ T. The approximation is sufficiently 
accurate as only a narrow energy level near the Fermi surface is responsible for heat conduction of strongly degenerate 
particles, m is related to the Debye mass (mo) through the relation m 2 = m 2 D /2. Hence, 7 can be expressed as, 

(33) 



f e2 1 








[12^ log 
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Thermal electron conductivity now takes the following form, 

&rlog 



12t 2 



1 4 ' 

^ 7T A . 



a 2 T~ 2 



7T° 



VC(3) 




(34) 

It is known that the specific heat capacity is c„ = /x 2 T/3 at the LO. For the thermal relaxation time we now have, 
from the Eq.([T]), 



9 



'0 



15 



(35) 



In the above equation the dominant first term in the second bracket is from the magnetic sector, the second term 
is from the longitudinal transverse cross term and the third one is from the electric one. This should be contrasted 
with the earlier results what have been reported in [f| @. In [j| @ the authors have mentioned that in the low 
energy transfer region (u> << mjj, T « tud) the thermal conductivity becomes independent of temperature. This 
happens if only the LO term in A in the Eg. ([33]) is considered, but the coefficients of the other NLO two terms make 
them comparable with the first one. Further modification of the result takes place when medium modified dispersion 
relation is included in the phase space factor of r Kee . Inclusion of 7 changes the temperature dependence significantly. 
In this context we can comment here that 7 has to be included in the expression of T Kee since at very low temperature 
a 2 log(T) term becomes large in comparison to one. We also observe here that the Fermi liquid description breaks 
down with the inclusion of the magnetic interaction. 



B. Higher orders in the thermal relaxation time 



In this subsection we evaluate the higher order correction terms in the low temperature thermal relaxation time. 
This has its origin in the inclusion of the NLO terms in the specific heat capacity and 7. A convenient starting point 
for this would be to extend the calculation of 7 in the Ea. ([3"Tj) which relates itself to the electron self-energy (£). In 
case of the quasiparticle momenta close to the Fermi momentum the £ is dominated by the soft photon exchange. 
Beyond the leading order at low temperature, it is given by the following equation Q , 



£(e,fc) 



e 2 m 



12tt 2 7 



loE 



V 7T£ / 
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6144 - 256tt 2 
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20 



24-7TTO 
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A , 775 — (sgn(e) - V3i) 



864tt 6 



45tt 7 /3 \ m ) 
(sgn(e) + V3i) 

0.928 m\ i7rsgn(e) 



lo 



O 



(36) 



The dominant logarithmic term in the fermion self-energy comes from the transverse sector and it gives rise to 
logarithmic singularity when — > /1 or in other words for excitations near the Fermi surface. We approximate the 
above expression in the low temperature limit as (e^ — fi) ^ T, as has been done in the last subsection. With the 
above quasiparticle self-energy we have, 

dek V duj J 



1 



127T 2 

(1 + 13), 



log 



4 

77a 



2 2/3 e 2 A 2/3 40 x 2 l/3 e 2 A 4/3 



9V3tt 7 /3 



27\/37r 11 /3 



(37) 
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where, j3 = 7 + 7 . 7 gives us the NLO NFL terms, 

, ^ 2 2 / 3 e 2 A 2 / 3 _ 40 x 2 1 / 3 e 2 A 4 / 3 
7 ~ 9V3?r 7 / 3 27V37r n /3 

The final expression for the electron thermal conductivity now becomes, 

1 



(38) 



|% 2 T-(l + /3 ) [A 2 C(3) + &f-X*<>< (| ) 



—A 3 
15 



(39) 



Unlike Fermi-liquid result where, K ee varies inversely with T, here the temperature dependence is non-analytical, 
anomalous in nature reminiscent of many other recent studies involving ultradegenerate plasma [H, [l2l . IbH . fl5l | . We 
show here that K ee involves fractional power in (T/m£>) coming from the medium modified phase space factor. 

The other quantity which we require for the estimation of relaxation time is the specific heat. It has already been 
derived in the context of degenerate quark matter in 0, [13]. For the degenerate electron gas it can be written as, 
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-40 



3(i 



2 2/3 r (i) c m TO 3 



ttA 

3 

3 J "*-D A 5/3 + 56(] 



2 1/3 r(f)C(f)^ A7/3 



27y/3nV 3 81V3W 11 / 3 
With the Eqs. (|3T))) and (|4U|) the relaxation time for thermal conduction is found to be, 



(40) 
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8lV37r n /3 



_D-)J/3 



(41) 



It is evident that the thermal relaxation time upto NLO terms contains some anomalous fractional powers origi- 
nated from the transverse interaction. This in turn changes the temperature dependence of r Kcc non-trivially. The 
appearance of the non-analytic terms in the Eqs.([39|). (j40|) and (|41|) has a common origin as explained earlier. 



III. RESULTS AND DISCUSSION 



In this section an estimation of the electron thermal conductivity and the electron thermal relaxation time with 
the temperature has been presented. In the left panel of the Fig.([T]) we note that the inclusion of both the HDL 
propagator and (3 decreases the value of K ee . It shows strong deviation from the Fermi liquid result n ee oc l/T. In 
the right panel it has been shown that j3 enhances thermal conductivity. This has serious implication on the total 
electron conductivity K e , since, the contribution of n ee has now increased in comparison with K ej ;. In [(| the authors 
have shown that magnetic interaction decreases n ee which in turn increases the electron-electron collision frequency. 
Thus to the total electron thermal conductivity electron-electron scattering dominates over electron-ion scattering. 
The phase space correction due to the medium modification of the electron dispersion relation partially compensates 
for such an enhancement. 

With the thermal conductivity and specific heat capacity we further plot the thermal relaxation time with the 
temperature. In Fig. ^ we have shown how the thermal relaxation time changes from the FL result with the inclusion 
of the medium modifications. In the right panel it has been shown that the inclusion of {3 enhances the thermal 
relaxation time. We have taken the relevant region of temperature and density from the reference 0. According to 
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FIG. 1: Temperature dependence of the electron thermal conductivity. The left panel shows a comparison between the Fermi 
liquid and the non-Fermi-liquid NLO result where /3 7^ 0. The right panel shows the enhancement of K ee with the inclusion of 



Q in the density region higher than p = 10 6 gem -1 and temperature more than 10 8 K Landau damping becomes 
important since in this density region non-relativistic degenerate electron gas becomes relativistic. This region is 
relevant for our plots. We have taken a — .3 and Dcbye mass 344 KeV, this satisfies the condition of T << mr> where 
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FIG. 2: Temperature dependence of the thermal relaxation time. The left panel shows a comparison between the Fermi liquid 
and the non-Fermi liquid NLO result when /3 7^ 0. The right panel shows the enhancement of r Kei! with the inclusion of /3. 



IV. SUMMARY 



In this paper, we calculate the thermal relaxation time in degenerate electron gas in the domain where the relativistic 
effects become important. It has been shown that with the inclusion of the magnetic interaction, which is relevant 
in the above mentioned domain, r Kce shows departure from the FL behavior. It is known that for the normal FL i.e 
with only Coulomb interaction r Kee behaves as 1/T 2 . In the relativistic domain, we on the other hand, find this to be 
proportional to the inverse temperature i.e t Kcc oc 1/T. This has been attributed to the absence of screening in the 
transverse sector. Wc further expose how the in medium modifications of the electron dispersion characteristic affects 
the heat conduction from the neutron star crust to the core. Our calculation actually modifies the phase space or the 
Fcrmionic density of states as revealed in the text leading to an enhancement of conductivity or relaxation time. For 
the thermal relaxation time, a closed form analytical expression has been derived by making perturbative expansion 
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in (T/mr>) retaining terms beyond LO. The appearance of the fractional powers in these results are interesting which 
is reminiscent of what one obtains in the calculation of the fermionic self-energy at high density. Numerically these 
corrections have also been found to be important in the present context. Particularly these corrections become large 
in the domain of small frequency i.e when uj << mn- 
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Appendix 



In this appendix we provide the complete details of the derivation of the energy integral in the phase space factor 
discussed in the main text. We substitute — (e p — fi)/T = rj p , — (e^ — fi)/T = rj^ and —lu/T = z to obtain, 

I ^/-( 1± /(^))( 1± /(— ))de P de k =T*J o — - - ^ - ^ — - - ^ - ^ d Vp d m . (42) 



The lower integration limit can be send to infinity without introducing much error. Now, following [20tl21| the general 
way to calculate the integral (for a function f = if, where s is an integer) as follows, 



where, 



2C 2 , = — J ^^^(l + sX-ire-^dr? 
(-1) 



t |s] is the integral part of and Co = g. Hence, one can write, 



f 00 z 

T ~, r^r-n ^ d V P = T- . (45) 

J (e-'fr + 1) (e^+ 2 ) + 1) p 1 - e" z v ; 

Finally, wc obtain the result of the phase space energy integration quoted in the Eq. (|31|) in the main text, 

de p de k f e J ek (l±f (€p+u) )(l±f {ek _ u) )=T 2 - Y (46) 
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